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^ ' In the dynamical triangulation model of 4D euclidean quantum 

gravity we measure two-point functions of the scalar curvature as a 
function of the geodesic distance. To get the correlations it turns out 
that we need to subtract a squared one-point function which, although 
<D \ this seems paradoxical, depends on the distance. At the transition and 

^ ' in the elongated phase we observe a power law behaviour, while in 

the crumpled phase we cannot find a simple function to describe it. 
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1 Introduction 



In the dynamical triangulation model of four dimensional euclidean quantum 
gravity the path integral over metrics on a certain manifold is defined by a 
weighted sum over all ways to glue four-simplices together at the faces |], ||]. 
This idea was first formulated in [§], using hypercubes instead of simplices. 
The partition function of the model at some fixed volume N is 

Z(N,k 2 )= J2 exp( K2 iV 2 ). (1) 

T(Ni=N) 

The sum is over all ways to glue iV four-simplices together, such that the resulting 
complex has some fixed topology which is usually (as well as in this article) 
taken to be S A . The iVj are the number of z-simplices in this complex. /t 2 is 
a coupling constant which is proportional to the inverse of the bare Newton 
constant: n 2 oc Gq . 

It turns out that the model has two phases. For low k 2 the system is in 
a crumpled phase, where the average number of simplices around a vertex is 
large and the average distance between two simplices is small. In this phase the 
volume within a distance r appears to increase exponentially with r, a behaviour 
like that of a space with constant negative curvature. At high k 2 the system is 
in an elongated phase and resembles a branched polymer. As is the case with a 
branched polymer, the (large scale) internal fractal dimension is 2. The transition 
between the two phases occurs at a critical value k\ which depends somewhat on 
N. This transition appears to be a continuous one, making a continuum limit 
possible |4], H, ||. At the transition, the space behaves in several respects like the 
four dimensional sphere 0. 



2 Curvature and volume 

In the Regge discretization of general relativity, all the simplices are pieces of 
flat space. The curvature is concentrated on the subsimplices of codimension 
two, in our case the triangles. On these triangles it is proportional to a two- 
dimensional delta function. From the definition of curvature as the rotation of 
a parallel transported vector, one can find the integrated curvature over a small 
region V t (A) around such a triangle 

f R^dx = 2A A 6 A , (2) 

JV e (A) 

where A A is the area of the triangle and 5 a is the deficit angle around the triangle 
(see e.g. ||). The deficit angle around a triangle is the angle which is missing 
from 2tt 

5a = 2tt- J2 °d (3) 

dG{S(A)} 
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where {S'(A)} are the simplices around the triangle and 9d is the angle between 
those two faces of the simplex that border the triangle. The angle S A can be 
negative. 

In dynamical triangulation, all the simplices have the same size and shape 
and the deficit angle is a simple function of the number ka of simplices around 
the triangle. Then expression (Q) reduces to 

/ R JT) dx = 2V 2 (2tt - 9n A ) , (4) 

JV e (A) 

where 9 is the angle between two faces of a simplex, which for D dimensions 
equals arccos(l/D), and V 2 = A A is the now constant area of a two simplex. 

For each triangle we can define a local four volume that belongs to the triangle 
by assigning that part of each adjoining simplex to it which is closer to the triangle 
than to any other. For equal simplices, this just results in V4/IO per adjoining 
simplex with V4 the volume of a four simplex. In other words this local volume 
is 

V A = ^/gdx = — (5) 
Jn(A) 10 

where 0(A) is the region of space associated to that triangle. It is not clear what 
Va would mean in the continuum limit. We define it here mainly to compare our 
results with other work on simplicial quantum gravity. 

If we view the delta function curvature as the average of a constant curvature 
over the region 0(A), this constant curvature would be equal to 

20V 2 2tt - 9n A 
Ra = t~t • (6) 

V4 TlA 

Because neither a constant term nor a constant factor is important for the 
behaviour of correlation functions, we will in the rest of this paper for convenience 
use the definitions 

Ra = n A \ (7) 
V A = n A . (8) 



3 Two-point functions 

One of the interesting aspects of the dynamical triangulation model one can 
investigate is the behaviour of two-point functions of local observables. Because 
we are looking at observables defined on the triangles we consider correlations 
between the triangles, at a fixed distance d which is also defined in terms of 
triangles. Such a correlation function of an observable 0(x) will be denoted by 
(00)(d). 
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We define the geodesic distance between two triangles as the smallest number 
of steps between neighbouring triangles needed to get from one to the other. For 
this purpose, we define two triangles to be neighbours if they are subsimplices 
of the same four-simplex and share an edge. Other definitions of neighbour 
are conceivable. One such a definition would be to define two triangles to be 
neighbours if they share an edge, irrespective of whether they are in the same 
simplex. The one we use has the advantage that it is quite narrow and therefore 
results in larger distances. 

The idea behind our correlation functions is as follows. For each configuration 
generated according to the ensemble ([!]), we take a random pair (x, y) of triangles 
at distance d, where x and y denote the triangles. For this pair we calculate 
the observable 0{x)0(y). Then we go to the next configuration and repeat the 
process. Finally we take the average over all such pairs. If no such pair exists for 
a particular configuration, the configuration is discarded. 

Obviously, this method would be very inefficient in practice. We improve the 
statistics by using the average value of 0(x)0(y) over all pairs (x, y) at distance 
d in each configuration. Therefore, we calculate for each configuration 



Y J 0(x)0{y)5 d ^ y)td 

x,y 



d(x,y),, 



(9) 



x,y 

Taking the average over configurations 



J2Aex V (K 2 N 2 (T)) 

{A) = ( N (TV ' (10) 

2^exp(K 2 iV 2 (T)) 

r 



results in the correlation function 



(00)(d) = (^— V (11) 



y)4 



Other definitions are conceivable. One can take random pairs of triangles from 
the collection of all configurations. Configurations with relatively many pairs of 
triangles at distance d will then be counted more often. In formula, it results in 
the correlation function 

(j20(x)0(y)5 d{Xj y )jd \ 
(OO)'(d) = Af* ^—L. (12) 

( £ $d(x,y),d ) 
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Figure 1: The correlation function (RR)(d) — (R) 2 for various values of k<i- 



A few experiments did not show a qualitative difference in the behaviour of (|TT|) 
and (|T2|) at the distances considered below. However, for large distances where the 
finite size of the configurations comes into play, the difference becomes significant. 

A third possibility, which is natural in Regge calculus, treats Tt x clS db local 
volume element at x in a discrete approximation to a continuum integral over 
euclidean spacetime: 

(y2n x n y O(x)0(y)5 d{Xty)4 ) 
(00)"(d) = r '-. (13) 

In this paper we explore the form (|Tl~l) . We expect that the correlations con- 
structed from either (|TTD, (|l^) or (|l"3|) will behave identically for not too large 
distances. In a large system, compared to the distance under consideration, the 
sum over the triangles will introduce a self-averaging which probably makes the 
difference in averaging between ([l]]) and (|12D irrelevant. 

In figure [I] we have plotted the correlation function of the curvature, with the 
square of its expectation value subtracted. Most of the data in this paper are for 
a volume N = 16000 simplices. The values of k,2 correspond to a system in the 
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Figure 2: The correlation function (nn)(d) — (n) 2 for various values of k 2 - 



crumpled phase (k 2 = 0.8), near (but slightly below) the transition (re 2 = 1.22) 
and in the elongated phase (k 2 = 1.5). 

Configurations were recorded every 10000 sweeps, where a sweep is defined as 
a number of accepted moves equal to the number of simplices N. For k 2 = 0.8, 
1.22 and 1.5 we used 16, 51 and 21 configurations respectively. 

One thing is immediately striking, the correlation functions do not go to zero 
at long distances. To keep the short distance behaviour visible, the full range in 
the elongated phase has not been plotted, but we already see that also in this 
phase it crosses the zero axis and indeed this curve does eventually go to large 
(« 0.02) positive values. 

The local volume Va is proportional to the number of simplices around a 
triangle i. We see that in this model, this observable Va is essentially the same 
as the scalar curvature. At first sight, one would therefore expect them to have 
the same behaviour. If one is positively correlated, the other one would also be 
positively correlated. Figure |2| shows the correlation of n. We see that quite the 
opposite is true. With few exceptions, n is positively correlated where R = n^ 1 
is negatively correlated and vice versa. 

This behaviour is similar to that reported for the Regge calculus formulation 
of simplicial quantum gravity in |J. There it is also found that the curvature 
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correlations are positive and the volume correlations negative at large distances. 

This difference in behaviour can be explained intuitively as follows. Because 
triangles with large n have more neighbours, any random triangle will have a 
large chance to be close to a point with large n and a small change to be close 
to a point with small n. So whatever the value of n at the origin, the points 
nearby have large n and the points far away have small n. The average (nn) will 
then be large at small distances and small at large distances. Because large n 
means small R, the situation is reversed if we substitute R for n in this discussion, 
qualitatively explaining figs. [I] and |[ 

At first sight one might conclude from this explanation that a point with 
large n having many neighbours is just an artefact of the model. This is not 
true, however. Large n corresponds to large negative curvature and also in the 
continuum a point with large negative curvature has a larger neighbourhood. To 
be more precise, the volume of d- dimensional space within a radius r around a 
point with scalar curvature R equals 

V(R) = C d r d (l - ^^r 2 + 0(r 4 )). (14) 

4 Connected part 

The above reasoning leads us to the somewhat unusual concept of a correlation 
function which does not depend on some observable at the origin. We define such 
a correlation as 

(R)(d) = (^- ), (15) 

\ l^ d d{x,y),d I 

where x and y denote a triangle. In the more usual case of a quantum field theory 
on flat space this could never depend on the distance, but here it does. The reason 
is that we correlate functions of the geometry with the distance, which is itself a 
function of the geometry. 

Figure |3] shows this correlation function. No average has been subtracted. 
The behaviour of this one-point function turns out to be very similar to that of 
the curvature correlation in figure |l|. This correlation function again shows that 
any particular point has a large chance to be in the neighbourhood of a point 
with low curvature, which can be simply explained with the fact that points with 
low curvature have more neighbourhood. 

The same plot for n (not shown) shows the opposite behaviour. At small 
distances it is larger than average, while at large distances it is smaller than 
average. This is rather obvious, because where n is large, its inverse is small and 
vice versa. 

We can now investigate how much of the curvature correlation shown of fig- 
ure [T] is due to this effect. Figure [| compares the curvature correlation (RR)(d) 
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Figure 3: The curvature as a function of the distance (R)(d) for various values 
of k 2 . 



with the square of this one-point function. We see that, except at small dis- 
tances, the two are indistinguishable on this scale. In other words, we have not 
been measuring any curvature correlations. All we have measured are correla- 
tions between the curvature and the geodesic distance. Similarly, (nn)(d) and 
{n)(d) 2 are nearly equal. 

It is now easy to explain the difference in behaviour between the curvature and 
the volume correlations. Because they are almost equal to the square of (R)(d) 
and {n)(d) respectively, they behave just like them. And as we just mentioned it 
is easy to understand that these have opposite behaviours. 

The way to go now is to subtract the two things and see what real curvature 
correlations are left. This is similar to subtracting a disconnected diagram and 
keeping the connected part. We get the corrected correlation functions 



C R (d) = (RR)(d) - (R)(d) 2 
C v (d) = (nn)(d) - (n)(d) 2 



(16) 
(17) 



The results for the curvature are plotted in figure [| and those for the volume 
in figure [f| The errorbars were found by a jackknife method, each time leaving 
out one of the configurations in the calculation of Cn(d) and Cy{d). Now both 
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correlations behave almost exactly the same. Note the large difference in scale 
between these figures and figures || and ||. 

In the crumpled phase we were not able to fit Cn(d) to a simple function. 
This is probably due to the fact that we cannot reach very large distances in this 
phase. Near the transition however it is possible to fit the correlation function to 
a power law decay, at not too small distances. This is shown in figure [7[ In the 
region 9 < d < 18 it fits nicely to ad b with the result 

a = -0.5(2) (18) 
b = -4.0(2) (19) 
X 2 = 5 at 8 d.o.f. (20) 

This data was made at a volume of 32000 simplices, with K2 = 1.255. We used 
65 configurations, which were recorded every 5000 sweeps. A similar fit using Cy 
gives a compatible power, 

a = -5.7(1.6) x 10 2 (21) 
b = -4.30(12) (22) 
X 2 = 2.3 at 8 d.o.f. (23) 
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Figure 5: Corrected correlation function Cr{<£) at various values of K2- 

This result should be taken with caution, however. One would really like 
to have a good fit over a larger range. To get some idea of the typical ranges 
involved, we consider the number of triangles at a given distance d 

(N'(d)) = I ^L— \ (24) 

where N2 is the number of triangles of the configuration. The corresponding 
quantity with 'triangles' replaced by 'four-simplices' was studied more closely in 
[[7]. The value d m where N'(d) has its maximum, is an indication of the distance 
at which finite size effects might become important. At k>2 = 1-255, this d m is 
only 11, indicating that finite size effects may play a role in the measured power. 

The situation is even better far in the elongated phase. Here a power law fits 
well, as can be seen in figure ||. This fit was done to the points 3 < d < 15 and 
the parameters of this fit are for Cr 

a = -0.0038(1) (25) 
b = -2.56(3) (26) 
X 2 = 17 at 11 d.o.f. (27) 
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Figure 6: Corrected correlation function Cy{d) at various values of k 2 - 
and for Cy 

a = -2.2(1) (28) 
b = -2.57(2) (29) 
X 2 = 14 at 11 d.o.f. (30) 

The value of d with maximum number of triangles was 32, so in this case 
we are in a region of small distances compared to the system size. This data 
was made from 23 configurations of 32000 simplices. We have also fitted the 
connected correlation functions at other points far in the elongated phase and at 
16000 simplices. The power that emerged was within the errors equal to the one 
given above. 



5 Discussion 

We have investigated the behaviour of the curvature and volume correlation func- 
tions. It turned out that the naive correlation functions could be almost entirely 
described by a "disconnected part", which we therefore subtracted. The differ- 
ence turns out to behave according to a power law in the elongated phase and 
near the transition. This indicates the presence of massless excitations. 
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Figure 7: Power law fit to curvature correlation Cr(<£) near the phase transition 
at iV 4 = 32000 and k 2 = 1.255. 



An obvious question is what the continuum observables are corresponding to 
the correlation functions we have measured. Generally, in a scaling region, a 
lattice operator is equivalent to a combination of various continuum operators, 
weighted with powers of the lattice distance a according to their dimensions. In 
a previous paper |7j we found evidence for scaling in a surprisingly wide region 
around the transition value of Assuming that the theory can be described by 
a continuum metric tensor g^ v with corresponding curvature R, the continuum 
observable with lowest dimension corresponding to our lattice correlation function 
would be given by 



(I dx^Jg(x) J dy^/g(y) 5(d(x, y) - d)R(x)R(y)) 
(I dx^fgjx) J dy^gjy) 8(d(x, y) - d)) 

(I dx^Jg(x) J dy^fg{y)S(d(x, y) - d)R{x)) 



(J dxJg(x)JdyJg(y)5(d(x,y) - d)) 



(31) 



where d(x, y) is the geodesic distance between the points x and y for a given 
metric g^ v . The uniqueness of the lowest dimension correlation function in the 
continuum is in accordance with the fact that we found identical behavior for Cr 
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Figure 8: Power law fit to curvature correlation Cr(cI) in the elongated phase at 
k 2 = 1.5. 



and Cy, up to an overall factor. Of course, we do not know the effective action 
specifying the average in (pT|). It could be a combination of / dx^fgR and higher 
order R 2 terms. 

The fact that the connected correlation function is negative suggests that it 
registers fluctuations in the conformal mode of g^. 

In a previous paper pj we explored the possibility of a semiclassical region 
near the transition, in which the system behaves like a four-sphere for not too 
small or large distances. To this end, we defined a scale dependent effective 
curvature. For k 2 near the transition the following picture emerged. At small 
distances, this effective curvature is large, indicating a Planckian regime. At 
intermediate distances there seems to be a semiclassical regime, where the space 
behaves like a four-sphere. The fluctuations around this approximate S A might 
then correspond to gravitons. We consider it therefore encouraging that we find 
the power law behaviour. 

For the volumes in current use, the effective curvature shows that the semi- 
classical regime sets in at a distance roughly 0.6 of r m (cf. fig. 13 in f?]). Here, 
r m is the geodesic distance through the simplices where the number of simplices 
N'(r) has its maximum. We conjectured this fraction to go down at larger vol- 



13 



umes. Similarly, a little beyond 0.6 of d m = 11 turns out to be the distance where 
the curvature correlations start to behave like d~ 4 in fig. [7]. We like to think of 
this as a confirmation of the point of view sketched above. 

Two-point functions of curvature and volume have been studied in the Regge 
calculus formulation of simplicial quantum gravity in refs. |9|, [TO]] . In these studies 
there are only results in what is called the well-defined phase of the Regge calculus 
approach, which corresponds to our crumpled phase. This makes it hard to do 
more than the qualitative comparison which was done in section |3]. 

The curvature correlations have also been investigated in the continuum. In 



ref. |TT|] they were found to be of zero range in the tree approximation to Einstein 



gravity. To one loop order we may expect on dimensional grounds a behavior 
G 2 d~ 8 in flat space, and G 2 R d~ 4 for S 4 , where G is Newton's constant and R is 
the S 4 background curvature. 

In [[HJ a theory is developed for the conformal factor in four-dimensional 
quantum gravity and from this the curvature correlation is calculated. The con- 
formally invariant phase discussed in [O seems to correspond to the elongated 



phase in the dynamical triangulation model. Intuitively, this can be understood 
by visualizing large fluctuations in the conformal factor as generating many baby 
universes. Many baby universes is also a feature of the branched polymer like 



elongated phase of simplicial quantum gravity fll3fl . Furthermore, the conformally 
invariant phase would occur at very large distance scales. In [7j we argued that 
the elongated phase also describes scales which are large compared to a typi- 
cal physical curvature scale. In this conformally invariant phase a power law is 
predicted for the curvature correlations (see also [14|]). Unfortunately, a direct 



comparison with |12], |1J] is not possible because in the continuum the correla- 
tion function is defined as a function of the distance in a fixed fiducial metric, a 
quantity that is not yet defined in our model. Our result « 2.6 for the power in 
the elongated phase is quite different from the ~ 0.7, which corresponds to the 



analogue central charge Q ~ 8 suggested in [14 
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